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SOME  REMARKS  CONCERNING  RECENT  WORK  ON  ROTATING  TURBULENCE 

YE  ZHOUt,  CHARLES  G.  SPEZIALE^,  AND  ROBERT  RUBINSTEIN^ 


Abstract.  A  recent  paper  on  rotating  turbulence  by  Canuto  and  Dubovikov  is  examined  from  both 
an  historical  and  scientific  perspective.  It  is  first  shown  that  their  claim  of  finding  a  new  energy  spectrum 
scaling  is  inaccurate;  such  a  scaling  law  has  been  published  in  the  literature  by  other  authors  using  the  same 
physical  assumptions.  Canuto  and  Dubovikov  actually  only  offered  a  different  estimate  for  the  constant. 
Finally,  it  is  demonstrated  that  the  alternative  model  for  the  dissipation  rate  transport  equation  proposed 
by  Canuto  and  Dubovikov  does  not  have  the  desired  physical  features  in  rotating  isotropic  turbulence.  It  is 
physically  inconsistent  in  both  the  weak  and  strong  rotation  limits. 

Key  words,  turbulence  modeling,  rotating  turbulence,  rotation  modified  energy  spectrum 

Subject  classification.  Fluid  Mechanics 

1.  Introduction.  Recently,  Canuto  and  Dubovikov  [1-6]  published  six  papers  under  the  general  title 
“A  dynamical  model  of  turbulence” .  An  application  to  rotating  turbulence  by  Canuto  and  Dubovikov  [5] 
(hereafter,  denoted  as  CD)  was  offered  by  the  authors  as  an  example  that  the  extension  of  their  “dynamical 
model”  is  capable  of  “providing  a  coherent  explanation  of  a  wide  variety  of  apparently  disjoint  qualitative 
and  quantitative  LES  results”. 

We  feel  that  it  is  appropriate  to  address  the  three  major  issues  in  CD: 

(1)  It  is  rather  regrettable  that  CD  made  somewhat  inaccurate  statements  on  what  is  new  and  what 
has  been  published  in  the  literature  by  other  researchers  in  the  field. 

(2)  The  abstract  of  CD  stated  that  for  rotating  turbulence,  “the  spectrum  exhibits  a  new  form  E(k)  ~ 
(eO)1/2*;”2. 

(3)  In  the  introduction,  the  authors  further  stated  that  “we  find  that  the  spectrum  exhibits  a  new  form 
E(k)  ~  (eO)1/2^”2  ”. 

These  statements  are  rather  misleading  and  give  the  readers  the  erroneous  impression  that  the  E(k)  ~ 
k~ 2  scaling  of  the  energy  spectrum  of  rotating  turbulence  was  first  found  by  CD.  Hidden  in  the  middle  of 
the  paper,  the  authors  admitted  that  “In  Ref.  [22]  a  spectrum  of  the  [same]  type  ...  was  obtained .” 

Indeed,  the  rotation  modified  energy  spectrum  has  been  published  by  Zhou  [7]  in  the  Phys.  Fluids.  Ref. 
[22]  of  CD  is  based  two  fundamental  assumptions:  first,  the  energy  transfer  is  local,  and  second,  the  time 
scale  of  the  triple  velocity  correlation  in  the  energy  flux  function  is  of  the  the  order  of  I/O,.  Here  is  the 
rotation  rate  of  the  fluid.  These  two  assumptions  lead  directly  to  the  fc“2  scaling  of  the  rotation  modified 
energy  spectrum7 

(1)  E(k)  =  C'n(f2e)1/2fc-2, 
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where  e  is  the  dissipation  rate.  The  constant  Cq  is  estimated  as  Cq  =  1/A  —  1.22  —  1.87  for  the  typical 
range  of  the  Kolmogorov  constant7.  These  results  are  supported  by  recent  mathematical  analysis  [8]  and 
direct  numerical  simulations  [9]. 

2.  Analysis.  Now  we  turn  our  attention  to  CD,  and  review  the  procedure  by  which  they  derived  the 
m  spectrum  for  strongly  rotating  flows.  It  turns  out  that  the  main  physical  assumption  in  their 

“dynamical  model”  is  “that  the  transfer  of  energy  is  mainly  local  in  character Another  assumption  is 
that  “in  the  presence  of  rapid  rotation ,  the  triple  correlation  decays  on  a  time  scale  of  the  order  of 
Hence,  these  two  basic  assumptions  in  CD  are  exactly  the  same  as  in  previously  published  results7.  As  we 
have  shown  in  Ref.  [7],  these  two  assumptions  are  all  that  are  needed  for  deducing  an  E(k)  ~  fc~2  energy 
spectrum  for  strongly  rotating  turbulence. 

It  is  clear  that  the  statements  in  CD  and  those  cited  above  in  (2)  and  (3)  are  inaccurate.  CD  should  not 
have  made  the  claim  that  “we  find  that  the  spectrum  exhibits  a  new  form  E(k)  ~  (eQ,)l/2k~2”.  Instead, 
CD  only  estimated  a  proportionality  constant  with  a  different  value  (y45/8  to  be  specific.). 

CD  also  presented  a  critique  of  the  modeled  dissipation  rate  transport  equation 

(2)  e  =  -Cg  |  D  |  e 

in  Rubinstein  and  Zhou  [10].  The  energy  equation  for  isotropic  turbulence 

(3)  K  =  -e 

and  Eq.  (2)  implies  that  in  a  rotating  isotropic  turbulence,  the  turbulent  kinetic  energy  K  approaches  a 
constant  and  the  turbulent  dissipation  rate  e  approaches  zero: 

(4)  e(t)  =  e(  0)e-fit 

(5)  Ko(t)  =  Ko(0)-^[l-e-nf]. 

CD  argued  that  since  in  a  A:-2  inertial  range  cut  off  at  a  scale  fco,  the  relationship 

(6)  K  =  CnVeCikQ1, 

holds  and  the  integral  scale  k^1  must  increase  exponentially.  We  shall  now  address  the  issue  raised  by  CD. 

We  stress  that  in  rotating  turbulence,  both  a  k~2  and  a  Kolmogorov  fc-5/3  scaling  region  can  and  do 
co-exist.  The  continuous  transition  between  these  regions  was  treated  heuristically  by  Zhou7.  At  sufficiently 
low  Rossby  numbers  (strong  rotation),  the  energy  spectrum  E(k)  is  theoretically  expected  to  follow  the 
E(k)  ~  k~ 2  scaling  law  for  k  >  ko(t).  Here  ko(t)  denotes  the  lower  edge  of  the  inertial  range.  The  E(k) 
profile  shape  for  k  <  k0(t)  is,  in  general,  problem  dependent.  The  function  kQ(t)  is  a  property  of  the  solution. 
To  simplify  our  discussion,  we  use  the  “split  spectrum”  model 

{Cks  if  k  <  k0 

CnVtek-2  if  fc0  <  k  <  fen 

CKe2/sk-V3  if  k>kn. 

Here,  the  parameter  kn  =  (f^/e)1/2  appeared  in  Refs.  [11]  and  [7]  but  was  used  in  Ref.  [5]  without  reference. 
The  parameter  kn  separates  the  inertial  range  modified  by  rotation  (k  <  kn)  from  that  of  the  traditional 
Kolmogorov  inertial  scales  ( k  >  kn).  In  decaying  turbulence,  e  — >  0,  so  that  kn  — ►  oo.  We  therefore  consider 
the  case  when  the  Kolmogorov  inertial  range  does  not  appear  in  Eq.  (7).  The  existence  of  the  range  of 
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scales  k  <ko  suggests  an  alternative  to  the  spectral  dynamics  proposed  in  CD’s  argument  against  the  model 
Eq.  (2).  Namely,  suppose  that  in  the  limit  of  asymptotically  long  time,  &o  approaches  a  constant.  Then 
the  k~ 2  range  will  disappear,  since  e  approaches  zero.  Unlike  the  scenario  proposed  by  CD,  the  turbulence 
kinetic  energy  can  nevertheless  approach  a  constant  as  required  by  Eq.  (5),  because  energy  can  be  trapped 
in  the  far  infrared  region  of  scales  k  <  fc0,  where  it  will  undergo  purely  viscous  decay.  The  single-point 
model  proposed  in  Eq.  (2)  therefore  does  not  entail  indefinite  exponential  growth  of  the  integral  scale.  Since 
single-point  modeling  entails  a  drastic  loss  of  information,  it  is  not  surprising  that  a  given  single-point  model 
can  be  consistent  with  several  quite  distinct  models  of  the  underlying  spectral  dynamics. 

CD  propose  an  alternative  dissipation  rate  transport  equation.  It  is  also  essentially  that  proposed  by 
other  authors  and  actually  does  not  have  the  requisite  physical  properties.  There  is  a  problem  at  both 
weak  and  strong  rotation  rates  as  well  as  at  both  high  and  low  turbulence  Reynolds  numbers.  In  the  strong 
rotation  limit  of  isotropic  turbulence,  where  ft  — ►  oo,  the  dissipation  rate  transport  model  proposed  by  CD 
reduces  to  the  form 

(8)  e  =  -C€  2  — 

where  the  coefficient  Ce 2  is  approximately  equal  to  2.84.  This  gives  rise  to  a  power  law  decay  [12] 

(9)  K  ~  t~a 

where  the  exponent  a  is  approximately  equal  to  0.54.  The  turbulence  Reynolds  number  Rt  =  K2jvt  actually 
undergoes  a  weak  power  law  growth 

(10)  Rt  ~  t0 

where  the  exponent  (3  is  approximately  equal  to  0.46.  It  is  now  well  established  that  the  cascade  is  so 
disrupted  by  a  rapid  rotation  that,  at  high  turbulence  Reynolds  numbers,  the  turbulent  kinetic  energy 
remains  approximately  constant  and  e  — ►  0  (the  turbulence  undergoes  a  linearly  viscous  decay  since  the 
phase  coherence  needed  to  cascade  energy  from  the  large  to  the  small  scales  is  scrambled;  see  Mansour, 
Cambon  and  Speziale  [13]).  Hence,  there  is  no  question  that  a  power  law  decay  for  the  turbulent  kinetic 
energy  and  dissipation  rate  is  physically  incorrect  in  a  rotating  isotropic  turbulence  at  high  turbulence 
Reynolds  numbers.  Furthermore,  it  is  even  incorrect  for  low  turbulence  Reynolds  numbers.  It  has  been 
shown  by  DNS  that,  at  low  turbulence  Reynolds  numbers,  the  turbulence  Reynolds  number  peaks,  after 
undergoing  a  dramatic  increase,  and  then  decays  -  an  effect  that  cannot  be  represented  by  the  weak  power 
law  growth  given  in  (10).  This  is  a  precursor  to  the  high  Reynolds  number  result  of  e  — *  0  and,  thus, 
Rt  — ►  00.  It  results  from  the  imbalance  between  the  production  of  dissipation  by  vortex  stretching  and 
the  leading-order  part  of  the  destruction  of  dissipation  caused  by  the  rapid  rotation  which  scrambles  the 
turbulence  [13].  A  power  law  decay  cannot  properly  describe  this  effect.  More  sophisticated  models  -  that 
account  for  non-equilibrium  vortex  stretching  -  are  needed  [13]. 

As  far  as  the  weak  rotation  limit  is  concerned,  the  dissipation  rate  model  of  CD  reduces  to  the  asymptotic 
form 

(11)  i  =  -Ce2^-Ce3Ktf 

where  Ce 2  and  Ce 3  are  constants  that  assume  the  values  of  11/6  and  1/10,  respectively.  This  equation  can 
be  integrated  analytically  yielding  e/K  going  as  the  tan  ct  to  a  power  where  c  is  a  constant.  Thus,  the 
solution  is  oscillatory.  DNS  and  physical  experiments  have  not  indicated  that  there  are  oscillations  in  the 


3 


turbulent  kinetic  energy  and  dissipation  rate  for  this  problem.  So  this  model  is,  unquestionably,  incorrect. 
Equation  (11)  -  which  actually  forms  the  basis  of  the  Hanjalic  and  Launder  [14]  model  -  has  been  integrated 
numerically  by  a  Runge-Kutta  numerical  integration  technique  (see  Speziale  and  Gatski  [15]).  It  was  found 
that  -  for  an  initial  condition  of  QK0/e0  —  0.496  which  is  not  that  far  removed  from  the  region  of  small 
rotation  rates  -  the  turbulent  kinetic  energy  grew  during  the  first  few  eddy  turnover  times  and  the  solution 
became  unrealizable  through  the  development  of  negative  dissipation  rates  [15].  Realizability  is  violated  by 
(11)  since  e/K  =  0  is  not  a  fixed  point  in  that  equation  unlike  in  the  classical  dissipation  rate  equation; 
this  allows  negative  values  of  the  turbulent  dissipation  rate  to  develop  (e  =  0  is  no  longer  an  invariant  plane 
in  dynamical  systems  terms).  Hence,  it  is  clear  that  the  asymptotic  form  (11)  proposed  by  CD  for  weak 
rotation  rates  is  fraught  with  problems. 

Finally,  it  should  be  noted  that  the  results  of  CD  are  only  applicable  to  rotating  isotropic  turbulence 
where  it  still  has  problems  as  shown  above.  CD  gave  the  erroneous  impression  that  their  results  were  more 
general,  which  is  simply  not  accurate. 

3.  Conclusions.  We  have  examined  a  recent  paper  on  rotating  turbulence  by  Canuto  and  Dubovikov 
[Phys.  Fluids  9,  2132  (1997)]  from  both  an  historical  and  scientific  perspective.  We  have  shown  that  their 
claim  of  finding  a  new  energy  spectrum  scaling  is  inaccurate;  such  a  scaling  law  has  been  published  in 
the  literature  by  other  authors  using  the  same  physical  assumptions.  Canuto  and  Dubovikov  actually  only 
offered  a  different  estimate  for  the  constant.  Finally,  we  have  demonstrated  that  the  alternative  model  for  the 
dissipation  rate  transport  equation  proposed  by  Canuto  and  Dubovikov  does  not  have  the  desired  physical 
features  in  rotating  isotropic  turbulence.  It  is  physically  inconsistent  in  both  the  weak  and  strong  rotation 
limits 
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